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Abstract. We report time-resolved optical flux and circular polarization spectroscopy of the magnetic DA white dwarf 
HE 1045-0908 obtained with FORS1 at the ESO VLT. Considering published results, we estimate a likely rotational period of 
Prat — 2.7 h, but cannot exclude values as high as about 9 h. Our detailed Zeeman tomographic analysis reveals a field structure 
which is dominated by a quadrupole and contains additional dipole and octupole contributions, and which does not depend 
strongly on the assumed value of the period. A good fit to the Zeeman flux and polarization spectra is obtained if all field 
components are centred and inclinations of their magnetic axes with respect to each other are allowed for. The fit can be slightly 
improved if an offset from the centre of the star is included. The prevailing surface field strength is 16MG, but values between 
10 and ~ 75 MG do occur. We derive an effective photospheric temperature of HE 1045-0908 of T eS = 10 000 ± 1000 K. The 
tomographic code makes use of an extensive database of pre-computed Zeeman spectra (Paper I). 
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1. Introduction 

Until a few years ago, magnetism among white dwarfs 
had been considered a rare phenomenon. A fraction of 
~ 5 % of all known white dwarfs had been confirmed 
to be magnetic, with fi eld strengths covering the rang e 
from -SOkG-lOOOMG 1 JWickra masing he & Ferrariol EoOO). 
Presently, the low-field tail of the known field strength dis- 
tribution is established by four objects in the kilogauss range 
(B a 2-4 k G), which is the current dete c tion limit for 8-m class 
telesc opes JFabrika & Valvavii] 1 19991 lAznar Cuadrado et"al] 

2004) . Recent studies suggest a much higher fractional in- 
cidence of magnetic white dwarfs (MWDs) of at least 10% 
for objects with surface fields exceeding 2 MG , and prob ably 
even more if low-fie ld objects are included (Lie bert et alJ2 003: 
ISchmidt et ai]l2003t and references therein). While a high in- 
cidence is mainly found for cool, old white dwarfs, it is inter- 
esting to note that a high incidence of (weak) magnetic fields 
has also been detected in central stars of planet ary nebulae, 
which are the direct progenitors of white dwarfs (Jor dan et al. 

2005) , and in subdwarf B and O stars JO'Toole et alJl2005). 



Send offprint requests to: F. Euchner, e-mail: 
f euchner@astro . physik . uni-goettingen . de 

* Based on observations collected at the European Southern 
Observatory, Paranal, Chile, under programme IDs 63.P-0003(A) and 
64.P-0150(C). 

1 1 MG = 10 s Gauss = 100 Tesla 



There is strong evidence that the high-field magnetic white 
dwarfs have evolved from main-sequence Ap and Bp stars. 
Low- and intermediate-field objects are thought to originate ei- 
ther from late A stars that fall just above the mass limit be- 
low w hich fossil fields a re destroyed in the pre-main-sequence 
phase jTout et al.l E004). or from main sequence stars of still 
later spectral type. 

Since there is no known mechanism to generate very sttong 
magnetic fields in white dwarf interiors, the fields are be- 
lieved to be fossil remnant s of previous evolutionary stages 
(Brai thwaite & Spruitl |2Q04). Modelling of the field evolution 
showed that the characteristic time for Ohmic decay of the 
lowest poloidal multipole components is long compared with 
the white dwarf evolutionary timescale (Wen dell et aljfl987t 
Cummin 3 l2002h . Higher-order modes do not necessarily de- 
cay faster, however, since they may be enhanced by nonlin- 
ear coupling by the Hall effect if internal toroidal fields are 
present ( Musl imov et alJ 19951) . This is consistent with the find- 
ing of significant deviations from pure dipole configurations 
dBurleigh et all 1 9991: Klaxted et alJl200ri IReimers et al.ll2004 
lEuchner et al.l2005b . 

The present Zeeman tomographic analysis of phase- 
resolved circular spectropolarimetry of the white dwarf 
HE 1045-0908 provides further evidence for strongly non- 
dipolar fields. In the first paper of this series, we have demon- 
strated the ability of our code to derive the field configuration 
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Fig. 1. Left panel: Flux spectra of HE 1045-0908 from June 1999 (1-5) and December 1999 (6-7). Spectra 2-7 have been shifted 
upwards by two flux units each with an additional unit between spectra 5 and 6. Right panel: Combined flux and circular polar- 
ization spectra of HE 1045-0908 from the June/December 1999 observations, covering approximately one-half of the rotation 
cycle. These spectra, which have been collected into three phase bins, have been used as input spectra for the Zeeman tomo- 
graphic procedure. For clarity, the uppermost two curves have been shifted upwards by one and two flux units, respectively. The 
quoted phases refer to case (i) with P rot = 2.7 h. 



of rotating MWDs using phase-resolve d flux and circular po- 
larization spectra (Euch ner etail 2002. henceforth referred to 
as Paper I). In this and follow-up papers, we apply our code to 
individual objects. 



HE 1045-0908 was discovered in t he Hamburg/ESO ob - 
jective prism survey for bright quasars JWisotzki et alj fl996). 
Subsequent optical spectroscopy at the ESO 3.6-m telescope 
revealed a rich spectrum of Zeeman-split Balmer absorption 
lines and confirmed the object as a magnetic DA white dwarf 
tReimers et alll993) . By fitting theoretical Zeeman spectra for 
a centred dipole with a trial-and-error method, the best match 
was found by these authors for 7W _ 9200 K ? flj,, = 31MG, 



2. Observations 

We obtained rotational-phase resolved circular spectropo- 
larimetry for the magnetic DA white dwarf HE 1045-0908 
with FORS1 at the ESO VLT UTl/Antu in June and 
December 1999. The dates and times of the observations as 
well as the number of exposures and exposure times are given 
in Tabled The spectrograph was equipped with a thinned, anti- 
reflection coated 2048x2048-pixel Tektronix TK-2048EB4-1 
CCD detector. For all observations, the GRIS_300V+10 grism 
with order separation filter GG 375 covering the wavelength 
range ~ 3850-7500 A was used with a slit width of 1" yield- 



and a nearly equator-on view. Schmidt et al 
quently obtained a sequence of five flux and circular polar- 
ization spectra of HE 1045-0908 over a duration of 1 h. The 
shape of the flux spectra in their observation sequence changes 
monotonically from almost vanishing to strong Zeeman fea- 
tures, whereas the variation in circular polarization is less pro- 
nounced. They estimated that the 1-h interval represented ei- 
ther one-quarter or one-half of a complete rotation cycle with a 
probable rotational period of P ro , 2-4 h. 



(2001) subse- Table 1. Dates and times for the spectropolarimetric observa- 



tions of HE 1045-0908 obtained at the ESO VLT (f exp : expo 
sure time, n: number of exposures). 



Object 


Date 


UT 




fexp (min) 


n 


HE 1045-0908 


1999/06/09 


22:55- 


-00:21 


20 


4 




1999/06/10 


00:23- 


-00:37 


14 


1 




1999/12/06 


08:27- 


-08:44 


8 


2 
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ing a FWHM spectral resolution of 13 A at 5500 A. We were 
able to reach a signal-to-noise ratio S /N 100 per resolution 
bin for the individual flux spectra. The instrument was operated 
in spectropolarimetric (PMOS) mode. The polarization optics 
consists of a Wollaston prism for beam separation and two su- 
perachromatic phase retarder plate mosaics. Since both plates 
cannot be used simultaneously, only the circular polarization 
has been recorded using the quarter wave plate. Spectra of the 
target star and comparison stars in the field have been obtained 
simultaneously by using the multi-object spectroscopy mode 
of FORS 1 . This allows us to derive individual correction func- 
tions for the atmospheric absorption losses in the target spectra 
and to check for remnant instrumental polarization. 

2.1. Data reduction 

The observational data have been reduced according to stan- 
dard procedures (bias, fiat field, night sky subtraction, wave- 
length calibration, atmospheric extinction, flux calibration) us- 
ing the context MOS of the ESO MIDAS package. In order 
to eliminate observational biases caused by Stokes parameter 
crosstalk, the wavelength-dependent degree of circular polar- 
ization V/I has been computed from two consecutive exposures 
recorded with the quarter wave retarder plate rotated by +45° 
according to 



V 

7 



< f°-r 
\f°+r 



r-r 
f°+r 



(i) 



where f denotes the ordinary and f e the extraordinary 
bea m (see the F ORS User Manual for additional information, 
IJehin et albooll) . 

Since there were noticeable seeing variations during the ob- 
serving run, we applied a correction for the flux loss due to the 
finite slit width of 1", using the measured FWHM of the object 
spectrum at 5575 A to estimate the effect of seeing and assum- 
ing a Gaussian intensity distribution across the slit. 

2.2. Data analysis 

In June 1999, we obtained a sequence of five exposures of 
HE 1045-0908 covering a time interval of 1.7 h, terminated by 
bad weather. This run yielded two independent circular polar- 
ization spectra (Eq.0. In December 1999, another "snapshot" 
of two additional exposures was secured, yielding another po- 
larization spectrum. All flux spectra are shown in Fig. [2 (left 
panel). The temporal change in the five June 1999 spectr a is 
similar to that seen in the data set of ISchmidt et alJ ll200ll) . In 
our data, the Zeeman features a re most promin ent at the begin- 
ning of the observations, in the ISchmidt et al J data at the end. 

The features at the beginning and at the e nd of our run resemble 

I • 1 

those at the end and the beginning of the Schmidt et al. run, re- 
spectively, i. e. the variation o f the Zeeman features in our run 
is reversed with respect to the Schmid tet alJ data. Our isolated 
observation in December 1999 also fits into this pattern. We es- 
timate that the phase of the strongest Zeeman features occurs 
between our spectra 1 and 2 in Fig. [2 Spectrum 5 approxi- 
mately corresponds to the phase with the weakest Zeeman fea- 
tures. The implied rotational period is P rot 2.7 h if the com- 




5000 



Fig. 2. Zeeman maximum ((/> = 0) flux fx (in units of 
10~ 16 erg crrT 2 s _1 A" 1 ) and circular polarization V/I of 
HE 1045-0908, plotted along with the theoretically pre- 
dicted field-dependent transition wavelengths (A-B curves) for 
Balmer absorption lines. The top panel also shows a theoretical 
circular polarization spectrum for a model atmosphere perme- 
ated by a homogeneous field of B — 16MG (shifted upwards 
by 0.07 units, with the horizontal dashed grey line denoting the 
zero polarization level). In the bottom panel, filled circles de- 
note unambiguous identifications of transitions. 



bined Schmid ~et alJ and present observations cover a full ro- 
tational period. While the substantial variation of the Zeeman 
features suggests this might be true, there is unfortunately no 
proof of such a connection of the phase intervals covered in 
the separate o bservations. Alternatively, it is possible that the 
Schmidt et al. and our data do not cover a full rotational period 
and P rot > 2.7 h. We refer to the former as case (i) and to the 
latter as case (ii). In case (ii), there are several possibilities for 
the phase intervals covered by our and the Schmidt et al. data 
which we discuss below. The large variation in the strength of 
the Zeeman features cannot arise in too short a phase interval, 
however, and we find that for periods in excess of 9 h an accept- 
ab le solution is n o longer obtained. We follow the suggestion 
of Schmi dTet alJ of a period in the 2-4 h range and adopt 2.7 h 
as the preferred period, but report on the consequences of as- 
suming a longer period below. 

In preparation of the analysis, we note that the flux spectra 
of December 1999 (spectra 6 and 7 in Fig. [3 are very simi- 



4 



F. Euchner et al.: Zeeman tomography of magnetic white dwarfs. II. HE 1045-0908 



lar in shape to spectrum 5 and provide an additional indepen- 
dent circular polarization spectrum which connects in phase to 
the June 1999 run. We collect spectra 1/2, 3/4, and 6/7 into 
three flux and circular polarization rotational phase bins. For 
case (i) with P mt 2.7 h, these bins are approximately cen- 
tred at rotational phases <p - 0.0, 0.25, and 0.5, where </> = 
refers to the phase of maximum strength of the Zeeman fea- 
tures in our June 1 999 observatio n. The case (i) concatena- 
tion of our and the ISchmidt et al.l data requires that the flux 
spectrum at <p — 0.75 resembles that at <p — 0.25. As a repre- 
sentative case (ii), we consider twice the rotational period and 
tentatively assign the three spectra to <f> — 0.0, 0.125, and 0.25. 
Our observations now cover a phase interval of only A</> = 0.25. 
Furthermore, we adopt P I0t 7.5 h or even 1 1 .3 h and assign 
the spectra to (p = 0.0, 0.09, and 0.18 or 0.0, 0.06, and 0.12 
with A(f> = 0.18 or 0.12, respectively. 

In Fig.[T](right panel), we present the three flux and circular 
polarization spectra which form the basis of our tomographic 
analysis. We refer to the mean of spectra 1 and 2 as "Zeeman 
maximum" {<p = 0.0) and to the mean of spectra 6 and 7 as 
"Zeeman minimum", which corresponds to <p — 0.5 for case (i) 
and to the smaller values given above for case (ii). 



3. Qualitative analysis of the magnetic field 
geometry 

The rich Zeeman spectra of HE 1045-0908 allow us to ob- 
tain insight into the magnetic field geometry already by 
the simple means of comparing the spectra with the ex- 
pected field-dependent wavelengths of the hydrogen transi- 
tions /l^fi^jjiejiceforthre^^ dForster et alJ 
ll984URosner et alJl984tlWunner et alJl985h . Fig.fjshows the 
wavelength range around HB of the Zeeman maximum spec- 
trum ((/> - 0.0) along with the A-B curves. Several transitions 
that can be immediately identified are marked by filled grey 
circles. This holds for the <x and n components in the flux spec- 
trum and the <x components in the polarization spectrum, while 
the circular polarization of the n component vanishes, indicat- 
ing a small viewing angle iff between the magnetic field direc- 
tion and the line of sight. The distribution of field strengths is 
sharply concentrated at ~ 16 MG, as demonstrated also by the 
fair agreement between the observed and the model polariza- 
tion spectrum shown at the very top of Fig. [2] This model spec- 
trum is calculated for a single value B = 16 MG, iff = 29°, and 
Teff = 10000K. Hence, the field over the visible hemisphere at 
this phase is ~ 16 MG and directed towards us. 

The identification of Zeeman transitions is not as easily 
possible in the Zeeman minimum spectrum (<p = 0.5), and we 
do not show the corresponding attempt of a quick analysis. 
Four transitions can definitely be identified, however, from al- 
most stationary parts of the A-B curves in Ha and HB cr + , 
and the corresponding field strengths span a range from 20 
to 60 MG. This simple analysis proves already that the field 
strength over the stellar surface varies by about a factor four, 
excluding simple field configurations like a centred or a mod- 
erately offset dipole. 
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Fig. 3. Observed circular polarization spectra for phases 
<p = 0.0 and 0.5 (grey curves) and best-fit synthetic spectra 
(black curves) for different parametrizations of the magnetic 
field geometry. From top to bottom: centred dipole (D ctrd); 
centred quadrupole (Q ctrd); off-centred dipole (D offs); off- 
centred quadrupole (Q offs); centred, non-aligned combination 
of dipole, quadrupole, and octupole (DQO na ctrd); off-centred, 
non-aligned combination of dipole, quadrupole, and octupole 
(DQO na offs). All curves except for the bottom one have been 
shifted vertically by suitable amounts in V/I, with the horizon- 
tal dashed lines indicating the respective levels of zero polar- 
ization. The quoted phases refer to case (i) with P lot = 2.7 h. 



4. Zeeman tomography of the magnetic field 

Theoretical wavelength-dependent Stokes / and V profiles of 
magnetized white dwarf atmospheres can be computed by solv- 
ing the radiative transfer equations for given B, iff, T e s, log g, 
and the direction cosine ft = cos §, where § denotes the an- 
gle between the normal to the surface and the line of sight. A 
synthetic spectrum for a given magnetic topology can be de- 
scribed by a superposition of model spectra computed for dif- 
ferent parameter values. Our three-dimensional grid of 46 800 / 
and V model spectra covers 400 B values ( 1-400 MG, in 1 MG 
steps), nine iff values (equidistant in cos iff), and 13 temperatures 
(8000-50 000 K) for fixed logg = 8 and fx = 1 (Paper I). This 
database allows fast computations of synthetic spectra for any 
given magnetic field configuration without the need to solve 
the radiative transfer equations each time. Limb darkening is 
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Table 2. Best-fit magnetic parameters for the different parametrizations of the magnetic field shown in Fig. [21 The uncertainties 
in the last digit are denoted by the values in brackets. A short description of the individual models is followed by the formal value 
of the reduced^ 2 (see text). All models assume a rotational period of 2.7 h, except for model (7), which has been computed for 
5.4h. 
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(1) D ctrd ( centred dipole, ^L, = 121.9) 

57(3) -32(1) 14(l) rC 30(1) - - - 

(2) Q ctrd (centred quadrupole, x]^ = 104.9) 

11(2) - - - -36(1) 14(1) 18(1) 

(3) D offs (off-centred dipole, fa = 29.7) 

20(3) -27(3) 41 (5) 29(2) - - - - - 0.08(1) 0.02(1) 0.27(1) 

(4) Q offs (off-centred quadrupole, xLa = 27.2) 

18(3) - - - re -49(4) 20(1) 20(2) - - - -0.13(1) 0.01(1) 0.06(1) 

(5) DQO na ctrd (non-aligned, centred combination of dipole, quadrupole, and octupole, xLa = 26.8) 
20(3) -12(2) 39(4) 17(2) -45(4) 36(3) 34(6) -19(1) 60(5) 27(3) 

(6) DQO na offs (non-aligned, off-centred combination of dipole, quadrupole, and octupole, x^, = 24.5) 

17(3) -16(2) 71 (5) 344(2) -36(3) 21 (3) 138(6) -18(2) 70(9) 115(7) 0.07(1) -0.08(1) 0.31 (2) 

(7) DQO na offs (non-aligned, off-centred combination of dipole, quadrupole, and octupole, P I0t = 5.4h, xLa = 26.3) 

32(4) -13(1) 31 (4) 1(1) -37(3) 24(3) 209(4) -25(3) 69(12) 245(10) -0.08(1) 0.07(1) 0.27(2) 



accounted for in an approximate way by the linear interpola- 
tion 

h(p)/h,n=i =a + b(x. (2) 

For the sake of simplicity, the temperature- and wavelength- 
dependencies of a and b have been neglected (see the discus- 
sion in Paper I). Fitting the measured absolute flux distribu- 
tion of HE 1045-0908 with model spectra computed using the 
full radiative transfer method, we find an effective temperature 
of T e s = 10000 + 1000K. For this temperature, model spectra 
computed as a function of yt suggest a = 0.53 and b - 0.47. 
In the subsequent tomographic analysis, these values of T e s, a, 
and b were employed and kept constant. 

Our Zeeman tomographic code requires an appropriate 
parametrization of the magnetic field, such that for every loca- 
tion r on the stellar surface the magnetic field vector B(r, a) can 
be computed depending on a parameter vector a = (a\ , . . . , Om) 
of M free parameters describing the field geometry. Best-fit pa- 
rameters are determined by minimizing a penalty function as 
a measure for the misfit between model and observation. For 
that purpose, we employed the C programming language li- 
brary evoC by Trint & Utecht (1994) 2 which implements an 
evolutionary minimization strategy. The penalty function we 
used is the classical reduced^ 2 

1 " (fj - sj(a)f 
= ^ h erf (3) 

with the input data pixels fj, the model data pixels Sj, and 
the standard deviations cr,. We used 1321 pixels per phase 
for the individual flux and polarization spectra each, yielding 
N = 7926 pixels in total (A = 3900-7200 A, AA = 2.5 A). All 

2 ftp://biobio.bionik.tu-berlin.de/pub/software/evoC/ 



phases have been equally weighted, and flux and circular po- 
larization have also been given equal weight. In order to esti- 
mate the statistical noise, a Savitzky-Golay filter with a width 
of nine pixels (corresponding to 20 A) has been applied to the 
observed spectra. Subsequently, the standard deviations crj en- 
tering Eq. Q have been computed from the differences be- 
tween the filtered and original spectra for wavelength intervals 
of 250 A. The standard form of^ 2 ed was used as a suitable rel- 
ative goodness-of-fit measure, but the unavoidable systematic 
differences between the observed and theoretical spectra pre- 
vent that anything near^ 2 ed 1 can be achieved. 

In order to avoid that such systematic differences influence 
the analysis of the narrower Zeeman structures, we adjusted 
the model flux spectra to the observed spectra at a number 
of wavelengths outside obvious Zeeman features. This proce- 
dure improves the fit of the flux spectra to the data but does 
not affect the polarization spectra. The wavelengths in ques- 
tion are marked by ticks at the top of Fig. ^ (right panel). Due 
to the finite exposure times, a model spectrum corresponding 
to a given observed spectrum should in principle be computed 
from several model spectra covering the respective phase in- 
terval. Our code is able to account for this "phase smearing" 
effect, but the need to compute the additional spectra slows 
down the minimization procedure so much that we decided 
against this approach. After having obtained best-fit parame- 
ters, we computed spectra including the effect and found no 
significant differences. The statistical errors of the best-fit pa- 
rameters have been computed using the method described in 
Zhan g et alJJl986l) . 

Details of the radiative transfer calculations, the synthetisa- 
tion of model spectra, the geometry adopted for the description 
of the magnetic field configuration, and the fitting strategy are 
given in Paper I. 
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4.1. Field parametrization 

In theoretical terms, a magnetic field of very general shape 
(with the constraints that it is curl-free and generated only in 
the stellar interior) can be described by expanding the scalar 
magnetic potential in spherical harmonics depending on the 
indice s / and m for the d egree and order of the expansion 
jGaufj|ll838l lLangellfl987l) . For each (/, m)-combination, two 
free parameters gf and hf are assigned, while only one pa- 
rameter, g9, describes the zonal components with m = 0. Thus, 
the number of free parameters for an expansion up to degree I 
is 1(1 + 2) = 15 (24, 35) for / = 3 (4, 5), increasing rapidly with 
the maximum degree /. Another property of multipole expan- 
sions is the dependence of the degree I required for an adequate 
description on the choice of the reference axis. Consider, e.g., a 
non-aligned superposition of a dipole and a quadrupole, which 
can be described exactly by four parameters (g®, g\, h\, g®) if 
the reference axis coincides with the axis of symmetry of the 
quadrupole. If the reference axis points in a different direction, 
a finite I allows only an approximate description. 

In order to ensure stable and fast convergence of our op- 
timization scheme, it is necessary to minimize the number of 
free parameters. We adopt, therefore, a hybrid model which 
implements a superposition of zonal (m = 0) harmonics only, 
disregarding the other tesseral components with m + 0. We al- 
low for arbitrary tilt angles of the zonal components and also 
for off-centre shifts. With this configuration, it is possible to 
describe fairly complex geometries with fewer parameters than 
in the truncated multipole expansion which includes all tesseral 
components (see Paper I for a detailed description). 

4.2. Results 

4.2.1. Case (i): P m = 2.7 h 

In an attempt to find the best-fitting field geometry for 
HE 1045-0908, we considered a sequence of parametrizations 
with increasing complexity. In Fig. [3] we compare the observed 
circular polarization spectra at <p = 0.0 and <p - 0.5 with such a 
sequence of model spectra. The best-fit parameters and the cor- 
responding xLa values are listed in Table |2] For the two sim- 
plest configurations (centred dipole and centred quadrupole) 
we obtained no satisfactory fit to the observations. This can 
be easily explained by the range of field strengths generated 
by these configurations, which is too large for <p = 0.0 and 
too small for <p = 0.5. It is interesting to note that the cen- 
tred dipole, which obviously provides the least adequate de- 
scription, is the only configuration that yields an inclination of 
i > 20°, whereas for all other configurations an inclination of 
i a 10°-20° is obtained. 

If an appropriate offset from the stellar centre is intro- 
duced for the dipole and quadrupole configurations, the pos- 
sible range of surface field strengths increases and an extended 
region with a nearly constant field strength of 16MG can be 
generated. Simultaneously, on the opposite stellar hemisphere 
a smaller high-field region with a steeper field gradient is cre- 
ated. As expected, the off-centred dipole and quadrupole mod- 
els match the observations better than the centred configura- 



tions, with the quadrupole model fitting better than the dipole. 
In general, however, these simple models are unable to produce 
adequate fits to all details at all phases simultaneously. 

The next steps in complexity of the field configuration are 
represented by the superposition of dipole, quadrupole, and 
octupole and the introduction of an off-centre shift. Dipole- 
quadrupole combinations were not successful and the inclusion 
of the octupole is essential. 

In a first attempt, we allowed the three individual compo- 
nents to be inclined with respect to each other, but not to be 
offset from the centre. The best fit with this field parametriza- 
tion matches the observed flux and polarization spectra well for 
all rotational phases (Fig. |4] top panel). The frequency distri- 
bution of field strengths extends from 10 to 70 MG and peaks 
at 16 MG. At (p — 0.0 the distribution of field strengths drops 
steeply towards lower and higher fields, while at <p = 0.5 the 
peak is less pronounced and the distribution is much broader, 
implying that fields up to 70 MG contribute significantly to the 
Zeeman spectra. The B-tfr diagram (Fig.|4] bottom left panel) 
shows that for the fields above 30 MG the sign of cos ip is re- 
versed compared with the most frequent field of 16 MG. The 
picture of the field geometry (Fig.0] bottom right panel) shows 
a high-field spot with B up to 70 MG superimposed on a low- 
field background of 10-20 MG. The field geometry on the visi- 
ble part of the stellar surface is quadrupole-like with field lines 
leading from the high-field pole to an "equatorial" band. The 
field strengths and orientations of the individual components 
are given in line (5) of Table |2 (see also Eq. 7 in Paper I). 
With 45+4 MG, the quadrupole is more than three times as 
strong as the dipole with 12 + 2MG. The three field compo- 
nents are more or less aligned, with the quadrupole inclined by 
only 11°, and the octupole by 22° with respect to the dipole. 
The slight inclinations of the individual components with re- 
spect to each other produce the required widening of the high- 
field spot. Basically, the field structure is that of an oblique ro- 
tator with an angle of ~40° between field and rotational axes. 
The lower right panel, labelled fi mx /RwD> indicates the max- 
imum radial distances to which field lines extend in units of 
the white dwarf radius. Distances beyond 10/?wd (black) may 
denote open field lines 3 . 

A slightly better fit is obtained if we allow for a common 
offset from the centre for all three components. As can be seen 
from Fig.[3](top panel), this additional freedom leads, in partic- 
ular, to improvements in the model circular polarization which 
we consider significant: (i) the steep rise at 4 170^-220 A for 
4> = 0.0; (ii) the dips at 4790^1870 A and at 5300 A for <f> = 0.5; 
and (iii) the continuum polarization in the 5200-6000 A range. 
The B-tfr diagram (Fig. [5J bottom left panel) shows an en- 
hanced frequency of field strengths around 60 MG for all three 
phases. For (f> = 0.5 we still see a pronounced decrease at 
70 MG, but, in contrast to the previous configuration, there is a 
small but significant contribution from fields of 70-76 MG and 



3 This information may not be relevant for HE 1045-0908 but is 
useful in studies of accreting white dwarfs, because it allows the iden- 
tification of regions that are potential foot points of field lines involved 
in channeled accretion. Accreting white dwarfs form a part of our pro- 
gramme and will be dealt with in forthcoming publications. 
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Fig. 4. Zeeman tomographic analysis of the magnetic field configuration of HE 1045-0908 using a centred, non-aligned com- 
bination of dipole, quadrupole, and octupole. Top: Observed (grey curves) and best-fit synthetic spectra (black curves). The 
uppermost two flux (circular polarization) spectra have been shifted for clarity by 2 and 4 (0.1 and 0.2) units in fx {VI I). The 
quoted phases refer to case (i) with P m = 2.7 h. Bottom left: B-ij/ diagram, bottom right: absolute value of the surface magnetic 
field, cosine of the angle ijj between the magnetic field direction and the line of sight, and maximum radial distance reached by 
field lines in units of the white dwarf radius (see text). 
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Fig. 5. Zeeman tomographic analysis of the magnetic field configuration of HE 1045-0908 using an off-centred, non-aligned 
combination of dipole, quadrupole, and octupole. Top: Observed (grey curves) and best-fit synthetic spectra (black curves). The 
uppermost two flux (circular polarization) spectra have been shifted for clarity by 2 and 4 (0.1 and 0.2) units in fx {VI I). The 
quoted phases refer to case (i) with P m = 2.7 h. Bottom left: B-ij/ diagram, bottom right: absolute value of the surface magnetic 
field, cosine of the angle ijj between the magnetic field direction and the line of sight, and maximum radial distance reached by 
field lines in units of the white dwarf radius (see text). 
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the same direction as the prevailing field of 16 MG. Figure [5] 
(bottom right panel) shows a field geometry that is similar from 
a global point of view, but reveals a more complex structure in 
the high-field region with two separate areas of opposite field 
direction. This general similarity is produced, however, by an 
arrangement of the field components which is entirely differ- 
ent from the previous model (see also the examples given in 
Paper I). The quadrupole is now inclined by 90° with respect to 
the dipole, and the octupole is not far from orthogonal to both. 
The shift is primarily in the direction of the dipole (z'-axis), 
and, hence, shifts the quadrupole perpendicularly to its axis. 
The field strengths quoted for this model in line (6) of Table |2 
refer to the unshifted components, and the final surface values 
can be calculated with some additional algebra. The shape of 
the region with field lines reaching beyond IQR WD (black) has 
changed from a circular spot to an arc. Field lines still end in 
a linearly extended region just visible at the stellar limb for 
= 0.0. 

All models in lines ( 1 ) through (6) of Table |2 have the 
property that the <p = 0.75 Zeeman spectra resemble those at 
<t> _ 0.25, as r equired by the proposed concatenat ion of the 
ISchmidt et al.l and our data. Hence, although the lSchmidt et alJ 
data have not been used in the fit, they are approximately re- 
produced by the models. 

The models of lines (5) and (6) of Table|2]with 9 and 12 free 
field parameters, respectively, fit better than the 17-parameter 
model of the full multip ole expansion up to I — 3 presented in a 
preliminary report feuchner et alJl2005l) . While this multipole 
expansion provides a better fit, e. g., to the 5200-6000 A con- 
tinuum polarization at <p — 0.0, it fails more seriously in other 
places. It seems that the choice of the arbitrarily oriented zonal 
components is more adequate for the case of HE 1045-0908. 
Models with about a dozen free magnetic field parameters rep- 
resent the present limit of our code at which a stable con- 
vergence to the global minimum in the x 1 landscape can be 
achieved. 

The general agreement between the observed and synthetic 
flux and polarization spectra has reached a high level which 
indicates that the theoretical spectra describe the underlying 
physics of magnetic atmospheres more or less correctly by now. 
The remaining differences can be traced back to a number of 
sources. On the theoretical side these are: (i) uncertainties in the 
absorption coefficients and approximations in the treatment of 
the line broadening, in particular, the treatment of Stark broad- 
ening in magnetic atmospheres; (ii) the finite resolution of the 
database, currently limited to 1 MG, which causes small wig- 
gles in the spectra. On the observational side these are: (iii) re- 
maining problems with the flux calibration, i.e. in the observa- 
tionally derived response functions, which we have attempted 
to correct for by re-normalizing the observed and model spec- 
tra relative to each other; (iv) small errors in the flat fielding 
procedure; and (v) uncertainties in the definition of the stan- 
dard deviations crj of spectral flux and polarization which enter 
Eq. and determine x] eA - 



4.2.2. Case (ii): P mt > 2.7 h 

We have repeated the analysis for rotational periods exceed- 
ing the preferred value of 2.7 h, in order to investigate whether 
the assumption of a longer period, which implies incomplete 
phase coverage, leads to a different field structure. The some- 
what surprising, but also fortunate result is that none of the 
investigated models yields a field structure which deviates sub- 
stantially from the one derived above. 

We replace the assumed value 0.50 of the phase interval A<p 
covered by our data by 0.25, 0.18, and 0.12, corresponding to 
rotational periods of 5.4 h, 7.5 h, and 1 1 .3 h. We consider first 
the case of P rot = 5.4 h. The important finding is that the as- 
sumption A<p = 0.25 does not imply the occurrence of a double 
wave of full period 5.4 h in the Zeeman features, but rather a 
field structure similar to case (i) seen at the larger inclination 
of i 32°. At least, this is true for our models which lack mul- 
tipole components higher than the octupole. As an example, 
we list in Table |2] line (7), the parameters for a non-aligned, 
off-centred dipole-quadrupole-octupole combination with an 
assumed period of 5.4 h, which should be compared with the 
model in li ne (6) for our preferred period of 2.7 h. For both 
models, the Schmi dt et alJ and our data are concatenated at the 
phase of Zeeman maximum, but in line (6) the combined data 
cover a full rotational period, and in line (7) only half a pe- 
riod. We conclude that the dominance of the quadrupole and 
octupole over the dipole is not affected by the different choice 
of the rotational period. 

The case (ii) model with A<p = 0.18 requires an inclination 
of i = 34°, whereas for A<p = 0.12 (with i = 53°) a satisfactory 
fit could no longer be obtained. This suggests a maximum value 
of the rotational period of about 9 h. It is not surprising that for 
decreasing A<p the necessary field variation between Zeeman 
maximum and minimum can only be produced by larger incli- 
nations which allow for a more rapid variation of the Zeeman 
features. 

5. Discussion 

In this study, we have fitted model Zeeman spectra to high- 
quality spectropolarimetric data of HE 1045-0908 using our 
Zeeman tomography code (Paper I), assuming a bona fide ro- 
tational period of about 2.7 h. We have achieved a good fit 
which reveals a dominant quadrupole component with addi- 
tional dipole and octupole contributions. HE 1045-0908 is the 
first white dwarf in which a quadrupole component has been 
detected so clearly. This result is found to be robust against the 
assumption of a longer rotational period with an upper limit at 
about 9 h. In our model, the orientations of the axes of dipole, 
quadrupole, and octupole have been treated as free parameters, 
and it turned out that this freedom is important in obtaining 
the best fit. This assumption deviates from a truncated multi- 
pole expansion with all m + components and is justified by 
its simplicity and ease of visualization. We are confident that 
we have reached a reliable reconstruction of the general field 
structure of HE 1045-0908. 

The most frequent photospheric field strength and direc- 
tion is represented by the maximum in the B—iff diagram at 
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16MG and positive cosip. This is also the field which appears 
most prominently in the Zeeman spectra, and a cursory anal- 
ysis would catalogue this star as "having a field strength of 
16 MG". Other sections of the star display field strengths up to 
~ 75 MG, however, which are less conspicuous in the observed 
spectra. Considering the complete information on the field dis- 
tribution, we find it difficult to assign either a "characteristic 
field strength" or a "polar field strength" to HE 1045-0908. 
More appropriate would be quotations like: (i) the most fre- 
quent field is 16 MG; (ii) the mean field over the visible sur- 
face averaged by the surface area is 34 MG; and (iii) the range 
of field strengths is 10-75 MG. Our general experience is, how- 
ever, that a quotation of type (iii) is model-dependent because 
models for some stars studied by us imply a high-field exten- 
sion in the B-ijj diagram, which covers only a small area near 
the limb of the visible surface and has little statistical signifi- 
cance. 

The derived field structure of HE 1045-0908 is primar- 
ily defined by the B-tft diagram rather than by the strengths 
and angles of the individual components. As pointed out in 
Paper I, different parameter combinations of the individual 
components could l ead to a similar B-tfr diagram. This is why 
iDonati et al.l {l994) refrained from specifying multipole com- 
ponents and su ggested to dir ectly optimize the B-iff diagram. 
The approach of lDonati et alJ does not guarantee, however, that 
the derived B-ip diagram corresponds to a physically possi- 
ble field. This potential trap is avoided in our approach, which 
has the additional advantage that we can specify the contribu- 
tions of individual multipole components. Furthermore, since 
we have gradually increased the level of complexity of our field 
parametrization starting from the elementary case of a centred 
dipole, we can be sure to have found the simplest configura- 
tion compatible with the observations. We cannot exclude ad- 
ditional small scale structure of the surface field, but suggest 
that such a structure cannot dominante HE 1045-0908 because 
it would destroy the remarkably high degree of circular polar- 
ization of up to ~ 10 %. 

Due to the small inclination (i = 17°) found for the best- 
fitting model geometry, a fraction of 35 % of the stellar surface 
is permanently hidden from view. In Paper I we have shown 
that this lack of information does not affect the accuracy of 
the derived field structure on the visible part of the surface. 
The field structures predicted by all our models for the hid- 
den part of the surface are reasonable and "well-behaved", i. e. 
there are no extreme field values or gradients. For instance, for 
the case (i) model of Fig. |5] and line (6) of Table [2] the range 
of field strenghts encountered on the visible fraction of the sur- 
face is 10-76 MG, while for the whole star it is 9-76 MG. Our 
simulations had shown that reconstruction artefacts can arise 
if the field parametrization involves more free parameters than 
needed to describe the field structure adequately (cf. Fig. 10 in 
Paper I). This does not happen in our stepwise approach with a 
limited number of parameters. 

It goes without saying that the final determination of the 
field structure of HE 1045-0908 would greatly benefit from a 
measurement of its rotational period and full phase coverage of 
the Zeeman spectropolarimetry. 



With r eff * 10 4 K, HE 1045-0908 is ~0.5Gyr old, less 
than the Ohm ic decay time s of the detected multipole compo- 
nents (Wendell et al. 1987; Cumming 2002). Hence, the strong 
quadrupole component could be a remnant from the main- 
sequence (or pre-main-sequence) evolution of the progenitor 
star. Alternatively, it could be pro duced by field evolution in the 
white dwarf stage as suggested by Muslimo v et alJ(ll995l) . Now 
that tomographic methods have reached a high accuracy thanks 
to advanced spectropolarimetric instruments at 8-m class tele- 
scopes, it would be interesting to follow up this question by a 
detailed field analysis of white dwarfs of different ages, supple- 
mented by more extensive theoretical calculations of the evo- 
lution of white dwarf magnetic fields over their cooling times. 
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